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ABSTRACT 

This  thesis  describes  briefly  the  scope  and  nature 
of  multi-echelon  inventory  systems  where  item  repair  is 
possible  at  the  various  levels.   A  model  is  created  which 
describes  many  commonly  encountered  real  world  systems. 
Sets  of  non-comparable  and  comparable  costs  are  specified 
and  a  cost-effectiveness  approach  to  the  solution  of  the 
model  is  outlined.   A  nonlinear  program  is  developed  where 
the  expected  sum  of  comparable  costs  is  minimized  subject 
to  a  given  level  of  fill  rate  (effectiveness)  provided  by 
the  lower  level  stocking  points  to  their  customers. 
Several  remarks  are  then  made  to  indicate  possible  solu- 
tion procedures  for  the  program  and  characteristics  and 
uses  for  the  cost-effectiveness  alternatives  developed. 
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I.   INTRODUCTION 

This  thesis  deals  with  a  complex  model  of  multi-echelon 
inventory  control  in  which  the  items  are  considered  to  be 
repairable.   Whenever  an  item  fails,  the  user  demands  a 
replacement  and  turns  in  a  nonfunctioning  unit  called  a 
carcass.   An  attempt  is  made  to  repair  the  carcass  and 
return  it  to  stock  for  future  issue.   Characteristics  of 
items  designated  as  repairable  are:   1)  such  items  are  phy- 
sically capable  of  repair,   2)  the  cost  of  repair  is  less 
than  procurement  cost,  and   3)  initial  procurement  costs 
are  large.   Complicating  this  scheme  is  the  fact  that  nearly 
all  such  models  are  of  the  multi-echelon  type.   That  is, 
several  levels  of  repair  and  stockage  interact  in  compli- 
cated ways  which  greatly  increase  the  difficulty  of  analy- 
sis.  Multi-echelon  systems  are  characterized  by  a  group  of 
stocking  activities  which  receive  customer  demands  and  in 
turn  demand  resupply  from  a  higher  level.   There  may  be  any 
number  of  levels  or  echelons;  however,  in  practice  only  tv/o 
or  three  levels  are  usually  encountered.   As  an  example  of 
a  multi-echelon  inventory  system  consider  a  large  retail 
store  chain  which  procures  and  stores  certain  types  of  goods 
nationally  and,  when  required,  supplies  regional  warehouses 
which  in  turn  supply  local  retail  outlets.   The  objectives 
of  the  study  of  multi-echelon  inventory  control  deal  to  a 
large  extent  with  the  development  of  those  po.l  Lcies  of 


ordering  at  all  levels  which  produce  some  specified  system 
performance,  such  as  minimum  expected  system  operating 
costs. 

In  addition  to  being  complex  multi-echelon  repairable 
item  inventory  systems  also  often  entail  huge  expenditures 
of  money  and  other  resources.   One  of  the  largest  multi- 
echelon  inventory  systems  dealing  with  repairable  items  is 
managed  by  the  U.  S.  Navy's  Aviation  Supply  Office  (ASO) . 
Repairables  represent  18.7%  of  the  items  controlled  by  ASO 
by  type,  for  a  total  of  66,000  separate  items;  but,  almost 
60%  by  dollar  value  for  an  investment  of  nearly  3.1  billion 
dollars.   The  Aviation  Supply  Office  spends  16  million 
dollars  per  month  for  repairs  at  the  lower  levels  of  repair 
alone,  that  is.  simple  or  relatively  inexpensive  repairs, 
and  28  million  dollars  monthly  for  procurement  of  such 
items.   Among  ASO  managed  repairable  items,  365,000  units 
were  reworked  in  fiscal  year  1971.   In  view  of  these  large 
investments  in  multi-echelon  repairable  item  systems  it  is 
not  difficult  to  understand  how  a  small  improvement  in  the 
way  in  which  these  systems  are  managed  can  result  in  sub- 
stantial savings. 

This  thesis  describes  a  model  which  approximates  the 
actual  working  structures  of  many  repairable  item,  multi- 
echelon  inventory  systems.   A  set  of  costs  and  measures  of 
effectiveness  is  defined  and  its  application  to  the  model  is 
discussed.   The  model  is  solved  for  minimum  expected  total 
system  cost  ordering  policy  subject  to  a  speci  ied  level  of 


effectiveness  or  performance.   The  solution  is  presented  in 
the  form  of  a  nonlinear  programming  problem.   Several 
remarks  are  then  made  in  an  attempt  to  outline  possible 
solution  procedures  for  the  program  and  applications  of  the 

model. 


II.   THE  MODEL 

The  lower  level  of  the  multi-echelon  repairable  item 
inventory  system  considered  in  this  thesis  will  consist  of 
two  bases  which  experience  demands  resulting  from  failures 
of  items  in  use.   When  a  demand  is  made  on  base  i,  a 
carcass  or  not-ready-for-issue  item  is  turned  in  at  that 
base's  repair  facility.   In  practice  it  is  the  case  that 
some  of  the  carcasses  which  are  turned  in  at  each  base 
cannot  be  repaired  at  the  base.   It  will  be  assumed  in  this 
thesis  that  an  item  will  be  repaired  at  base  i  with  proba- 
bility Pj:  .   If  repair  at  base  i  is  possible  the  repaired 
item  is  returned  to  stock  at  base  i  after  some  time,  rt., 
later.   It  is  assumed  that   rtj   is  a  random  variable  with 
some  general  distribution  and  has  mean  r^-±-       The  repair 
policy  at  base  i  may  then  best  be  described  as  an  infinite 
server  queue.   If  an  item  is  deemed  to  be  unrepairable  at  a 
base  it  is  forwarded  to  a  higher  level  repair  facility, 
depot  repair.   Define  the  quantity  of  items  on  hand  in 
stock  and  in  repair  plus  items  on  order  less  items  ordered 
by  customers  but  not  supplied  at  base  i  to  be  the  base  i 
inventory  position.   When  the  base  i  inventory  position 
reaches  a  level   r-j_,  a  request  is  made  to  the  depot  for  an 
additional  unit  which  arrives  a  time   t^   later  where   tj_ 
is  a  value  of  a  random  variable.   Since  the  base  i  inventory 
position  only  decreases  whenever  a  demand  is  received  and 


the  accompanying  carcass  is  not  repairable,  an  order  is 
placed  only  when  such  a  demand  occurs. 

A  single  depot  receives  demands  and  carcasses  from  the 
bases  when  repair  is  not  possible  at  base  level.   It  is 
clear  that  some  items  will  eventually  be  lost  since  their 
condition  will  have  deteriorated  beyond  repair.   For  this 
reason  it  is  assumed  that  carcasses  received  from  base  i 
may  be  repaired  with  probability   q.   and  if  repair  is  not 
possible  they  are  discarded.   Items  accepted  for  repair  are 
assumed  to  be  returned  to  depot  stock  after  a  constant 
repair  time,  rt.   Define  the  quantity  at  depot  both  on  hand 
and  in  repair  to  be  the  depot  net  inventory.   Due  to  the 
attrition  caused  by  items  being  discarded,  the  depot  net 
inventory  position  will  fluctuate.   Define  the  quantity  on 
hand  in  stock  and  in  repair  plus  items  placed  on  order  by 
the  depot  less  orders  placed  by  the  bases  to  the  depot  but 
remaining  unfilled  to  be  the  depot  inventory  position. 
When  the  depot  inventory  position  reaches  a  level   r,  an 
order  is  initiated  by  the  depot  for  a  quantity   Q   from  a 
source  of  manufacture  or  higher  level  of  supply  which 
arrives  after  a  constant  lead  time   t.   This  policy  is  simi- 
lar to  the  well  known  (r,Q)  model  of  single  stock  points 
which  has  many  practical  advantages  as  well  as  intuitively 
appealing  logic  and  simplicity.   Among  its  advantages  is 
the  ease  with  which  it  may  be  understood  and  applied  to 
nearly  any  inventory  requirement.   It  is  necessary  to  con- 
sider the  items  in  depot  and  base  repairs  wher  developing 
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an  inventory  ordering  policy.   In  this  thesis  this  is 
accomplished  by  including  the  quantity  in  repair  in  the 
inventory  position.   In  this  manner  it  is  possible  to 
insure  that  the  quantity  of  items  in  the  entire  system 
remains  bounded.   The  movement  of  items  within  the  model  is 
illustrated  in  Fig.  1.   A  listing  of  parameters  and  vari- 
ables used  in  the  description  and  solution  of  the  model  is 
included  in  Table  I. 

A.   MAJOR  ASSUMPTIONS 

Major  assumptions  about  the  model  required  in  its 

solution  are: 

Base  i  experiences  demands  which  may  be  closely 
approximated  with  a  Poisson  distribution  with 
rate   Lj_. 

The  depot  repair  facility  completes  repairs  on 
items  in  a  constant  specified  time  not  greater 
than  the  procurement  lead  time  t  (when  repair 
is  possible) . 

The  procurement  lead  time  for  the  depot  is  a 
constant  value   t. 

Bases  do  not  experience  crossed  orders,  i.e. 
orders  are  received  from  the  depot  in  the  same 
sequence  in  which  they  were  placed. 

The  depot  is  never  required  to  place  an  order 
to  a  higher  level  of  supply  when  a  previous 
order  has  not  been  received,   This  assumption 
does  not  usually  distort  the  model  since  for 
the  usual  case  where  ordering  costs  are  high 
and  holding  costs  relatively  low,  orders  are 
placed  infrequently. 

An  order  quantity  of  one  is  appropriate  to  a 
minimum  expected  cost  solution  for  base  i.   As 
argued  by  Feeney  and  Sherbrooke  £l]  ,  where 
holding  costs  at  the  lower  level  are  typically 
higher  than  at  the  higher  level  because  of 
such  factors  as  restricted  storage  capacity 


and  ordering  costs  are  much  lower  at  the  lower 
level  than  at  the  higher  level  since  they  only 
entail  administrative  and  handling  costs 
instead  of  for  example,  contract  negotiation 
costs;  a  policy  which  orders  small  amounts 
often  and  thus  maintains  low  stock  levels 
usually  compares  very  favorably  with  the  mini- 
mum expected  cost  policy.   This  is  especially 
true  where  integer  quantities  only  are  con- 
sidered. 


B.   PARAMETERS  AND  VARIABLES 

A  few  comments  in  explanation  are  appropriate  for  some 
of  the  parameters  and  variables  listed  in  Table  I.   I.  , 
*ir'   *   an<^   -'-r   are  inventory  holding  costs  for  a  base  i 
stock,  base  i  repair,  depot  stock  and  depot  repair  respec- 
tively.  These  figures  represent  the  costs  composed  of  such 
things  as  storage,  handling,  loss,  breakage,  etc.  for 
holding  each  item  of  inventory  for  a  given  length  of  time. 
It  is  assumed  that  the  holding  costs  are  expressed  in  units 
of  dollars  per  item-year.   Therefore,  the  instantaneous 
rate  of  incurring  carrying  charges  can  be  expressed  as   Ix 
where   x   is  the  on  hand  inventory  level.   Although  diffi- 
cult to  specify  in  practice  and  often  not  varying  linearly 
with  stock  on  hand,  as  assumed  here,  it  is  considered 
necessary  and  acceptable  to  introduce  holding  costs  in  this 
manner. 

The  order  cost  for  the  depot,  A,  is  taken  to  be  total 
cost  associated  with  placing,  following  up  and  receiving  an 
order.   Components  of  A  might  be  the  costs  of  letting  a 
contract  and  receiving  the  order  at  the  warehouse. 
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Li 

Pi 

rti 
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mean  rate  of  demand  for  the  item  experienced  at 
base  i 

probability  that  a  carcass  turned  in  at  base  i 
may  be  repaired  at  base  i 

expected  repair  time  for  items  repaired  at  base  i 
mean  procurement  lead  time  experienced  at  base  i 
base  i  reorder  level 

inventory  holding  cost  for  base  i  stock 
inventory  holding  cost  for  base  i  repair 

probability  that  an  item  forv/arded  to  the  depot 
?rom  base  i  repair  can  be  repaired  at  the  depot 

constant  repair  time  for  items  repaired  at  depot 

procurement  lead  time  experienced  at  depot  -  a 
constant 

r  depot  reorder  level 

q  depot  reorder  quantity 

inventory  holding  cost  for  the  depot  stock 
inventory  holding  cost  for  the  depot  repair 

A  order  cost  for  the  depot 

immediate  fill  -  major  components  of   d±   are 
shipping  time  and  time  required  for  processing 
of  papers,  packaging  for  the  orcier 


rt 
t 


I 
Ir 


Tab 


le  I.   Parameters  and  Variables 


12 


III.   SOLUTION  OF  THE  MODEL 

Solution  of  an  inventory  model  for  an  optimal  policy 
is  often  a  very  difficult  task.   First,  it  is  necessary  to 
determine  the  structure  of  the  problem  and  what  information 
is  available  and  desired.   Usually  the  analyst  proceeds  as 
though  the  consequences  or  dollar  costs  of  certain  states 
of  the  system  modeled  can  be  determined;  examples  are  the 
costs  of  holding  an  item  in  stock  without  issuing  it  for  a 
year  and  the  cost  of  placing  an  order  for  some  specified 
number  of  articles.   When  all  such  states  can  be  assigned  a 
cost  the  investigator  then  tries  to  determine  the  proba- 
bility of  the  occurrence  of  each  state  of  the  system  and 
then  the  expected  cost  of  maintaining  the  system  icr  a 
given  period  of  time.   It  is  necessary  that  this  expected 
cost  expression  be  a  function  of  variables  which  may  be 
determined  by  an  inventory  policy;  otherwise,  expected 
operating  costs  are  independent  of  inventory  policy  and 
there  is  no  problem.   Most  commonly,  these  policy  variables 
are  lead  times,  order  quantities  and  reorder  levels.   It  is 
then  necessary  to  solve  the  expected  cost  expression  for  a 
policy  which  yields  the  minimum  expected  cost.   Of  course, 
expected  costs  may  not  be  of  great  importance  and  it  might 
be  of  value  to  design  a  policy  which,  for  example,  guaran- 
tees with  a  certain  probability  that  all  demands  received 
can  be  filled  within  a  given  period  of  time.   The  analogy 
holds  in  this  case  also.   The  consequences  of  each  state 
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of  the  system  in  terms  of  the  probability  of  issue  within 
the  given  period  of  time  must  be  determined  and  the  proba- 
bility of  filling  an  order  within  the  given  period  of  time 
expressed  as  a  function  of  variables  which  the  inventory 
policy  controls.   As  in  the  case  with  expected  costs  it  is 
then  necessary  to  solve  this  problem  for  a  "best"  policy 
which  in  this  case  is  the  one  which  produces  the  desired 
probability  that  issues  are  made  within  the  given  time 
period.   It  is  important  to  realize  that  the  mathematical 
problem  which  results  from  the  type  of  formulation  discussed 
above  is  unconstrained;  that  is,  any  policy  which  produces 
the  desired  value  of  the  single  criteria  is  acceptable. 

A  problem  occasionally  arises,  however,  where  it  is  not 
possible  to  assign  costs  to  each  possible  state  which  are 
compatible  with  costs  associated  with  the  characteristics 
desired.   For  example,  while  an  investigator  may  attempt  to 
find  policies  which  minimize  total  expected  system  dollar 
costs  over  some  period  of  time  to  operate  a  system,  it  is 
not  always  possible  to  assign  costs  to  factors  such  as  not 
being  able  to  fill  a  customer  demand.   What,  for  instance, 
is  the  cost  in  dollars  of  not  having  an  item  in  stock  to 
repair  a  ship  damaged  at  sea  or  a  meteor  struck  space  craft? 
in  these  situations  the  best  that  the  analyst  can  usually 
hope  for  is  to  "pass  the  buck"  as  far  as  assigning  costs  of 
being  out  of  stock  is  concerned  to  the  authority  responsi- 
ble for  accomplishing  the  mission  involved,  i.e.  the  deci- 
sion maker.   In  this  formulation  the  analyst  ;  oduces  for 
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the  decision  maker  a  series  of  alternatives  of  expected 
costs  against  probability  of  being  able  to  fill  a  demand. 
The  alternatives  are  usually  arrived  at  by  solving  a  group 
of  problems  where  expected  value  of  the  sum  of  the  compar- 
able costs  is  minimized  subject  to  constraints  on  the  non- 
comparable  costs  or  consequences;  in  this  case  the  proba- 
bility of  not  being  able  to  fill  a  demand.   It  is  then  up 
to  the  decision  maker  to  select  a  policy  which  corresponds 
to  his  opinions  or  feelings  about  the  trade  off  between 
total  system  expected  cost  and  payoff,  or  in  the  parlance 
of  systems  analysis  he  must  specify  his  preference  ordering 
on  the  set  of  measures  of  cost  and  effectiveness.   If  the 
analyst  is  able  to  deduce  or  is  supplied  with  an  expression 
for  the  decision  maker's  trade  off  policy  between  cost  and 
effectiveness  in  this  situation,  he  may,  of  course,  apply 
this  policy  to  the  list  of  alternatives  he  has  developed 
and  determine  a  "best"  or  "optimal"  policy.   It  can  be 
shown  [Hadley  and  Whitin,  1963]  that  in  picking  a  policy  or 
formulating  a  preference  ordering  the  decision  maker  is  in 
effect  placing  a  dollar  value  cost  on  not  being  able  to 
fill  a  demand  and  in  this  sense  the  two  formulations 
discussed  above  may  be  considered  equivalent. 

In  this  thesis  a  formulation  similar  to  the  non- 
comparable  costs  example  discussed  above  is  used.   First  it 
is  necessary  to  specify  a  set  of  costs  to  be  considered;  or 
stated  differently,  it  is  necessary  to  determine  those 
characteristics  of  a  state  of  the  system  which  may  be 
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assigned  costs  which  are  comparable.   Other  characteristics 
which  are  important  but  which  are  not  included  above  are 
then  treated  as  constraining  factors  and  labeled  effective- 
ness criteria.   A  set  of  points  in  the  cost-effectiveness 
space  inferred  above  is  determined  by  the  solution  of  a 
program: 

Minimize  total  expected  system  costs 
(comparable)  subject  to  constraints  which 
require  that  each  measure  of  effectiveness 
(non-comparable  cost)  meets  or  exceeds 
some  stated  value 

for  many  different  levels  of  the  constraints.   Each  of  these 
points  is  simply  an  alternative  presented  to  the  decision 
maker  as  in  the  example  noted  above.   It  is  then  necessary 
to  apply  the  decision  maker's  preference  order,  if  avail- 
able., and  solve  the  related  problem  for  an  "optimal"  policy. 

This  relatively  elaborate  optimization  scheme  is 
developed  in  order  to  circumvent  the  very  difficult  task  of 
specifying  the  cost  of  not  being  able  to  fill  an  order 
immediately,  the  backorder  cost.   It  is  an  attempt  to  place 
the  introduction  of  this  consideration  into  the  system  model 
directly  through  the  decision  maker,  via  his  preference 
ordering.   Therefore  measures  of  effectiveness  must  neces- 
sarily address  the  manner  in  which  backorders  are  to  be 
treated  and  more  particularly  the  rate  of  backorders  experi- 
enced.  A  measure  of  effectiveness  is: 

B-L(ti,ri)  -  the  probability  that  base  i  is 

able  to  fill  a  demand  immediately  from  stock 
on  hand.   This  has  been  referred  to  by 
several  authors  as  the  fill  rate. 

There  is  one  such  measure  of  effectiveness  fo:  each  base. 
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The  total  expected  system  cost  to  be  minimized  is 
considered  to  consist  of  the  following  comparable  costs: 
total  expected  order  cost  for  the  depot  and  total  expected 
holding  costs  summed  over  bases  and  the  depot.   Although 
these  components  of  total  system  cost  are  also  difficult  to 
assess  in  certain  situations,  they  do  not  present  the  same 
possible  elements  of  "operational  imperative"  as  the  back- 
order  state  as  the  examples  of  the  damaged  ship  and  space 
craft  illustrate.   For  this  reason  and  because  of  the  diffi- 
culty of  arriving  at  an  optimal  policy  via  a  preference 
ordering  or  decision  maker's  judgment  when  the  set  of 
measures  of  effectiveness  is  large  (the  decision  maker  is 
offered  too  many  alternatives)  these  costs  are  assumed 
specified  and  included  in  the  objective  equation  of  total 
system  cost.   Mathematical  expressions  for  the  components 
of  the  program  are  developed  below. 

A.   CONSTRAINTS 

Let   B-^  (t^,  r^)   be  the  probability  that  base  i  is  able 
to  fill  a  given  demand  immediately  from  stock  on  hand. 
Define   uj_   to  be  the  number  of  items  in  repair  at  base  i 
plus  those  items  on  order  at  base  i,  v:   to  be  the  number 
of  items  in  repair  at  base  i  repair  and  w^   to  be  the 
number  of  items  on  order  at  base  i.   Then,  since  the  inven- 
tory position  of  each  base  is  a  constant   r^  +  1  under  the 
one  for  one  ordering  policy  used  by  the  bases: 
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Bi(t.,r-)    =   P (quantity    in  base    i   stock > 0) 


ri 


=  y       P(ui   =    J) 
j=0 


ri        J 

V      y        P(    v±    =    q,    w±    =    j    -   q) 
j=0      q=0 


ri 


=  y        y       Jp(v±    =    q)    p(Wi   =    j    -   q)J-   . 
j=0      q=0 

This  result  follows  from  a  finding  presented  by  Richards  (_3J 
which  shows  that  the  distributions  of  the  number  of  units 
in  repair  and  the  number  of  units  on   order  are  independent. 
A  summary  of  his  results  follows.   Item  failures  or  demands 
placed  upon  the  system  in  an  interval  of  time  (u,  u+s^ 
constitute  a  Poisson  process   <F(s),s  -  0 >  .   When  a  failure 
occurs  the  item  returned  to  base  i  is  designated  base 
repairable  with  probability   p^   and  is,  therefore,  not 
base  i  repairable  with  probability   l~Pi  •   it  i-s  assumed 
that  the  classification  of  failed  items  as  base  repairable 
or  not  base  repairable  is  independent  of  the  process  gener- 
ating failures  and  is  also  independent  among  failed  items. 
In  addition,  it  is  easily  shown  that  the  number  of  items 
designated  base  i  repairable  in  some  interval  of  time 
(u,u+s3,  Fi(s),  is  Poisson  with  mean  rate   LiP^   and  the 
number  of  items  which  are  not  base  repairable  in  that  length 
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of   time    is    Poisson  with   mean   rate      Li(l-pi).      Also,    the 
distribution   of   the   number    of    items   base    i   repairable    in   a 
length  of   time   given    that      n      total   demands   have  been  made 
in    the    same    length  of    time    is    simply  binomial  with  para- 
meters     n      and      p^.      With    these    facts    in  mind    it    is   only 
necessary   to   show  that   the   number   of   items  base    i   repair- 
able   received    in   a  period   of    time    is    independent  of   the 
number   of    items   which   are   not  base    i   repairable   received    in 
that    same   period   of   time.       If   this    is    true,    it    is   well   known 
that    the   number   of    items    in   each  classification   received    in 
different,    possibly   overlapping,    time    intervals    is    also 
independent.       Consider    the    following: 

p     |f2(u+s)     -   F2(u)    =   k,    Fx    (u+s)     -   Fx(u)    =   m^ 

=  p     Jf2(u+s)     -   F2(u)    =   k  J  F    (u+s)    -   F(u)    =   k+mV  ' 

p         F(u+s)     -   F(u)    =    k+m  \ 

Hin\        /  \  v  -L-js  .k+m 

me1       (Lis) 


k    ;       II    -  p±  1      p. 


(k  +   m)  ! 


k  m      -  PiLiS       -    (l-pi)Lis 

LU-Pi)    LisJ        (PjLjs)  e e 

k:    m! 

T/i  \T       "1  k      -(l-pj)L-iS  m      -p-jL-jS 

[_(l-pi)LisJ        e  ^i      i  (PiL±s)       e   ^i    ! 


kT 


m. 


p      )f2(u+s)     -    F2(u)    =   kip      Fx(u+s)     -    Fx(u)     =   m 
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Because  the  number  of  items  on  order  at  any  time  varies 
only  with  the  receipt  of  demands  accompanied  by  items  which 
are  not  base  repairable  and  the  procurement  lead  time,  the 
independence  of  the  distributions  of   Vj   and  Wj   is  esta- 
blished.  Because  of  an  important  theorem  by  Palm  [4J  it  is 
also  known  that  the  steady  state  probability  distribution 
of  the  number  of  items  in  repair  is  Poisson.  distributed  with 
parameter   rtiLiPi   an<^  the  number  of  items  on  order  is 
Poisson  distributed  with  parameter   t^L • (1-p • ) .   It  is 
interesting  to  observe  that  these  distributions  depend  on 
the  lead  times  only  through  their  mean  values.   Thus 


^  J-   f  ) 

i(ti,r1)=>"  )_        U^rtil^Pi)  ptJ-q.-t^d-Pi))^ 


j— u   q=u 

Where 


'. 


e 
p(x;A  )  =  ~ 


"X(A)X 


xl 


B.   OBJECTIVE  FUNCTION 

Let  K(r/Q/r1,r2)  be  the  total  expected  system  cost  for 
a  year.   Similarly  let   Ki(ri/ti)   and   Kd(r,Q)   be  the 
expected  costs  for  operating  base  i  and  the  depot  respec- 
tively.  Then   K(r/Q,r1,r2)  =  K^r^t^)  +  K2(r2,t2)  +  Kd(r7Q) 

1.   Development  of  the  Depot  Expected  Cost  Expression 


Let 


SD(r,Q)    =    E (number   of    items    in    stock   at   the   depot) 

RD  =   E (number   of    items    in   repair    at   the   depot) 

RSD(r,Q)=    E(inventory    status) 

OD(r,Q)    =    E(numbcr   of   orders    per   year) 

BD(r,Q)    =    E (number   of  backorders    at    the   depc   i) 
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where   E(X)   denotes  the  expected  value  of  the  random 
variable   X,  and  the  inventory  status  is  defined  to  be  the 
number  of  items  in  stock  and  in  repair  at  the  depot  less 
any  backordered  items. 

It  can  be  shown  that: 

Kd(r,Q)  =  (A)  00(1,0)  +  (Ir)RD  +  (I)SD(r,Q). 
Since  the  average  annual  demand  experienced  by  the  depot 
is  for   L-,  (1-q-, -p-, )  +  L^Cl-q^'Pp)  ~  L   items  the  expected 
number  of  orders  placed  by  the  depot  in  a  year  is   L/Q 
where   Q   is  the  size  of  the  order  placed.   It  may  also  be 
shown  using  Palm's  theorem  that  the  expected  number  of  items 
in  depot  repair,  RD,  is  simply   L(rt).   Since  the  inventory 
status  includes  the  number  of  items  on  hand  plus  those  in 
repair  less  backorders,  SD(r,Q)  -  RSD(r.O)  -  PD  +  BD(r,Q). 
In  addition  it  can  be  shown  that   ED(r,Q)  =  LP.   with 
probability  one  where   Pout-   is  t^e  probability  that  the 
number  of  items  on  hand  is  less  than  or  equal  to  zero. 
]?out   is  developed  in  the  next  section.   Therefore,  it  is 
only  necessary  to  determine  the  expected  value  of  the 
inventory  status,  RSD(r,Q).   This,  however,  is  exactly  the 
expression  developed  by  Hadley  and  Whitin  [_2J  for  the 
expected  on  hand  inventory  for  an  inventory  model  which  has 
no  repairable  feature  and  experiences  demands  which  are 
distributed  Poisson. 

Let  Z  denote  the  depot  inventory  position.  It  can  be 
shown  that  Z  may  assume  values  r+l,...,r+Q.  Z  does  not 
equal   r   for  any  finite  length  of  time  since  when   Z  =  r+1 
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and  a  demand  is  received  (recall  this  may  be  interpreted  as 
the  case  when  a  demand  is  received  v/ithout  a  repairable 
carcass)  an  order  is  placed  immediately  and   Z  =  r+Q.   It 
is  also  known  jGalliher,  Morse  and  Simon,  1959J  that  the 
inventory  position,  Z,  is  uniformly  distributed  on  the 
values   r+1, . . . . , r+Q,  i.e.   Z   takes  on  values   r+j   with 
probability  1/Q  where  j=l,2,...A. 

Let   X  be  the  inventory  status.   Consider  the  system 
at   u  -  t  where   t   is  the  procurement  lead  time.   Any 
order  placed  after   u  -  t  will  arrive  after  time   u.   Any 
order  placed  at  or  before   u  -  t  will  arrive  not  later 
than   u.   Then;  if   Z  =  r+j   at  time   u  -  t   the  probability 
that  there  are   x   units  on  hand  at   u   is  the  probability 
that   r  +  j  -  x  items  were  demanded  in  time   t   if   r  +  j  - 
x  -  0   and  zero  if   r+j  -  x  (  0.   Since  demand  is  Poisson 
distributed  with  parameter   Lt   the  probability  that 
r  +  j  -  x   demands  are  received  in  a  period   t   is 


-Lt    ,      ,  r+j-x 
e  (Lt) 

(r  +   j    -  x)  I 


if      r+j    -  x  )  0 


0  ifr+j-x<0 


Then 


Q        f  1 

(X  =  x)    =  V      jp(X=x  I    Z  =   r+j)    p(Z  =   r+j)J 


Q 

I 

j  =  l 
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Q 

=  §  Y~ p(r+:i  ■ x?  Lt)  if  _  °° <  x  <  r+i 
j=i 

Q 
=  —     V       P(r+j    -   x;    Lt)       if      r+1   <  x  <    r+Q. 

j=x-r 

Therefore: 

r+Q 

RSD(r,Q)    =      y       p(X  =   x)  . 
x=-oo 

It  then  follows  that: 

r+Q 
■SD(r,Q)  =   Y         p(X=x)  -  L(rt)  +  LPQut. 
x  =  -co 

2.   Development  of  Base  i  Expected  Cost  Expression 

By  an  argument  which  exactly  parallels  the  one 
above  for  the  depot  it  can  be  shown  that: 


Ki<ri.  ti)  -  Iir  Pi  L±  (rtj  +  I±1 


r  r-L+1 


/_  P(X±=X) 

L.  x=  -&Q 


-   p.  L.   (rt±)  +  L±      (1  -  P^t.P 


r-outj^ 


This  is  exactly  the  expression  for   K^(r,Q=l)   where  the 
term  corresponding  to  the  order  cost  is  omitted  since  it 
does  not  depend  on  the  variables   r.,  r   and   Q,  and  the 
demand  rate  and  procurement  lead  time  parameters  are  those 
appropriate  for  Base  i. 
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C.   DEVELOPMENT  OF  THE  EXPECTED  BASE  i  PROCUREMENT  LEAD 
TIME,  ti 

It  is  necessary  to  determine  an  expression  for   t^,  the 
expected  Base  i  procurement  lead  time,  which  appears  in  the 
expressions  for  KiCr^t^)  and   B^t-^,  r^  .   The  procurement 
lead  time  for  Base  i  can  be  decomposed  into  two  parts.   The 
first,  called  the  order  and  shipping  time,  is  that  time 
which  is  required  to  notify  the  depot  of  the  base's  need 
plus  the  actual  time  required  to  transport  the  unit  from 
the  depot  to  the  base  once  a  unit  becomes  available.   This 
time  will  be  assumed  to  have  expected  value   djj   and  is 
independent  of  the  inventory  situation  at  the  depot.   The 
second  part  is  the  amount  of  time  which  the  depot  requires 
to  obtain  a  unit  to  issue.   If   S   denotes  the  latter  time, 
then  it  is  clear  that   S   will  be  zero  if  and  only  if  the 
on  hand  inventory  at  the  depot  is  positive.   Thus   S   is  a 
random  variable  which  depends  on  the  on  hand  inventory  at 
depot,  the  items  in  repair  and  the  items  on  order  at  depot. 

The  expected  Base  i  procurement  lead  time  is,  therefore 


t.  =  d.  +  E(S) 
1     i     v 


Define  the  following  terms: 

W  =  the  random  number  of  items  in  depot  repair 


H  =  the  random  number  of  items  in  depot  supply 
ready  for  issue  less  any  backorders 


Y  =  W  +  H 
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Rs    =    the   number   of   receipts   by    the    depot    from 

procurement  and  repair  in  any  time  interval 
(u,u+sj  of  length  s 

Xs  =  the  number  of  receipts  by  the_  depot  from  repair 
in  any  time  interval  (u,u+s]  of  length  s 

Z   =  the  number  of  receipts  by  the  depot  from 

procurement  in  any  time  interval  (u,u+sj  of 
length  s. 


It  remains  to  determine  an  expression  for   E(S)   as  a 
function  of  the  program  variables   r   and   Q   to  completely 
specify  the  program  statement  sought.   For   S  -  0,  it  is 
true  that 

p(S  ^  s)  =       p(S  <  s  |  H  =  -h)  p(H  =  -h)  +  p(H  >  0) 


h=0 


where 


and 


p(H  =  -h)  =   ^   p(Y=x-h)  p(W  =  x) 
x=0 


'out  =   E   P  [  H  =  "h  1  ■ 

h=0 


As  developed  in  the  preceding  section  it  is  known  that 

Q 
p(Y=k)    =  ^      Y         P(r+j    -  k    ;    lit),  k    <    r+1 

3=1 

Q 

=  1     V  p(r+j    -   k;    Lt)  ,    r+1  £  k  <  r   +    Q 

j=k-r 

in    the    steady    state.       Therefore, 
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r+h     m  Q 


p(H  =  -h)  =  ^     )  q  Z   P  [  r+j  "  Cx-h);Lt]' 

x=0     ^   j=l 

p  f  x;  (rt)  L  ]  } 


r+h+Q 

/• 

Q 

1 

f* 

L 

x=r+h+ 1 

t 

j=x-h-r 

+  I      )h       I      P  Lr+j  "  Cx-h)lLt] 


p  [x;  (rt)  £]  I 


where   L  =  L^Q-j_  +  L2^2* 


Also  as  shown  in  the  preceding  section  the  random  variables 
X   and   Y   are  independent.   Therefore,  it  is  possible  to 
write 

p(H>   °)    =        )  p(H=h)    =      \  >  p(Y-h+x)    p(W=x) 

h=l  h=l        x=0 


r  r-h     f  q 

I    I  \h  £  p(r+j  - tx+h):  Lt)' 

h=l  x=0     <-         j=l 


r+Q  r+Q-h        ,  Q 

p(x;     (rt)L))        +       Y.  E  Q      L 

J  h=r+l      x=r+l-h  *»         j=: 


p(r+j    -(x+h);Lt)    p(x;     (rt)  L    )     f    . 


r+Q  r+Q-h         r         J2 

x-h-r 


J 
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It  remains,  therefore,  to  develop  an  expression  for   P(S  ^  s 

|  H  =  -h) 
Case  1;   Let   s  <  t.   (Recall  that   t  denotes  the  depot 
procurement  lead  time.)   Whenever   H  =  -h,   S   will  be  less 
than  or  equal  to   s   if  and  only  if   Rs   is  greater  than  h. 
Thus, 

p(S  i  s  |  H  =  -h)  =  p(Rs  >  h)  =  Y,  P(RS=J)  • 

j  =h+ 1 

Recall  that  a  necessary  consequence  of  the  assumption  that 
no  more  than  one  order  can  be  unfilled  at  any  time  is  that 
only  one  order  can  be  received  in  any  time  interval  of 
length   s  <  t.   Thus, 

p(Re=J)  =  p(Z-=0,  Xo=j)  +  p(Zs=Q,  X~  =  i-0). 

Continuing, 

p(Zs  =  i,  Xs  =  j)  =  p(Zs  =  i)  p(Xs  =  j) 

since  the  random  variables   Zg   and   Xs   are  independent. 

Proceeding, 

P(Z  =0)  =  P(no  order  arrives  in  an  interval  (u,u+sj  of 
length   s) 

=  P(no  items  are  on  order  at  time   u+s-t) 
r+Q 

I 

k=  r+1 


=  Y         P(Y  =  k) 
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r+Q 

=  y 

k=r+l 
p(Zs=Q)  =  1  -    P(ZS  =  0)  . 

Subcase  1;   Let   s  <  rt.   It  is  true  that: 

P(XS  =  j)  =         P(XS  =  j  |  W  =  x)  p(W=x). 
x=j 

Then 

P (X  = j   W=x)  =  P(x-j   items  are  in  repair  at  time  u+s 
each  of  which  was  in  repair  at  time   u) 

An  expression  for  the  steady  state  probability  that   x 
items  are  in  repair  each  of  which  has  been  in  repair   s 
units  of  time  can  be  developed  for  any  repair  time  distri- 
bution as  follows   J  Takacs  1962  j  .   Let 

a(s)  =  P(any  item  in  repair  in  the  steady  state  has 
been  in  repair  at  least   s   units  of  time) . 

Then  from  the  distribution  of  age  for  a  renewal  process  it 
is  known  that 


a(s)  =   jl/£    /  [l  -  F(t)]   dt 


where   T   is  the  mean  repair  time  and   F(t)  is  the 
probability  distribution  function  for  the  repair  time. 
Let   D   be  the  event   "y  units  are  in  repair   t  time   u+s 
each  of  which  was  in  repair  at  time   u"  and  lc  t   X  be  the 
number  of  items  in  repair  at  time   u.   Then 
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(PC 


P(Dy)   =    ^   P(DY  |  X=x)  P(X=x) 


X=j 


-   £   (y)  -<S)Y  d-a(s))X-y  <Ls>  x^±S 
a(s)y  (Ls)y  e"Ls   y  £ls  (l-a(s))]  *"* 


(x-y) : 


/v 


=   a(s)y  (Ls)y  e       _Ls(l-a(s)) 


A  -y   -Ls  a(s) 


(a(s)  Ls)1    e  "B  n^'         a 
— i =   p(y;Lsa(s)) 


Then, 


1  A 

P(X  =  j    W  =  x)   =   P(x-y;  Lsa(s)). 


For  the  case  where  the  repair  time  is  constant,  the  cumula- 
tive distribution  function  is  given  by: 

F(u)  =1      if   u  ^  rt 

=  0      if   0   ^  u  <  rt. 

Therefore   a(s)  =  1  -  s/rt. 

Subcase  2:   Let   rt  -  s  <  t. 

Since   rt  -  s   everything  in  repair  at  time   u  will  have 
been  repaired  by  time   u+s.   In  addition,  any  items  which 
enter  repair  in  the  interval  (u,u+s-rtj  will  also  have  been 
repaired.   Therefore,  the  conditional  probability  that  the 
number  of  items  received  from  repair  in  the  interval 
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(u,u+sj   is   j   given  that   x   items  are  in  repair  at  time 


u   is 


p  (Xs=j  |  W=x)  =  p(j-x;  (s-rt)  L  )    if   j  ^   x  . 

For   the   case    in  which      j   <  x. 

p(Xs=j  |  W=x)    =   0. 
Therefore,    for      0   -    s    <   rt: 


oo         s      r+Q 


p (S^s  |   H=-h)    = 


Q 

V 


j=h+l   (.     k=r+l 


E"  to       I         P(r+i-*;    Lt)] 
i=k-r 


oo 

Y_        P  [x-J*    Lsa(s)]      p    [x;  (rt)    Lj  + 


x=j 


r+Q 

\  ' 


Q 


.1.        \ 


1 


oo 


?,  .  T      ■  .  \ 


/  'q        £        Plrr:i-    k;    lj^) 


k=r+l  i=k-r 


x^j-Q 


p  [x-j+Q;    Lsa(s)J    p  [x;     (rt)    L  ] 


"\ 


^ 


And    for      rt   -    s  <  t: 


p(S^s   I  H= 


oo  f  r+Q  Q 

-w  =  Y.     <      L   \h    L     p(r+i~k; 

j=h+l      i      k=r+l    ■-       i=k-r 


Lt) 


p    [j-x;(S-rt)   £]     p    [x;  (rt)    L  _ 


x=0 
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r+Q        Q 
1  "   X.     Q    Z    P(r+i~k'  Lt) 


k=r+l 


i=k-r 


J-Q 

L 

x=0 


fj-Q-x;  (s-rt)  LJ  p  [x;  (rt)  L    J 


Case  2:   Let   s  ~  t. 

Because  it  is  assumed  that  not  more  than  one  order  is 
outstanding  at  any  one  time  it  follows  that: 

p  (S="s  |  H=-h)  =  1. 
With  the  expression  for   P(S-s)   thus  developed  it  is 

possible  to  calculate: 

t 

E(s)  =   i  p  (S  j>  s)ds  . 


D.   STATEMENT  OF  THE  PROGRAM 

The  minimum  expected  cost  solution  for  given  values  of 
b^   and  b2   where  h^      is  the  probability  that  an  order 
occurring  at  Base  i  can  be  filled  immediately  may  be  deter- 
mined by  solving  the  nonlinear  program: 


Minimize 


2 

T 

1=1 


Ki(ri'   V   +   Kd   (r'Q) 


subject  to 


Bi(t±/  r±)  £  bi#  i  =  1,  2 


where   r.,r   and   Q  may  only  assume  non-negative 


integer  values. 
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IV.   SOLUTION  OF  THE  PROGRAM 

This  thesis  does  not  develop  an  efficient  procedure  with 
which  the  program  specified  can  be  solved.   However,  several 
observations  are  offered  below  in  an  attempt  to  delineate 
the  scope  and  nature  of  the  difficulties  involved. 

A.  PROGRAM  TYPE 

It  should  be  noted  that  the  program  developed  is 
nonlinear  in  both  the  objective  function  and  the  constraints 
and  the  program  independent  variables  are   r,  ,  r~ ,  r  and  Q, 
each  of  which  may  assume  non-negative,  integer  values.   It 
is  not  apparent  that  the  program  is  convex,  and  it  is  almost 
certain  that  establishing  that  the  program  is  or  is  not 
convex  in  the  objective  equation  and/or  constraints  would 
involve  laborious  computations. 

B.  DEVELOPING  A  SOLUTION  PROCEDURE 

It  is  important  to  note  that  for  given  values  of   r   and 
Q,  r-^  and  r2   are  determinate,  i.e.  they  assume  the  smallest 
integer  values  such  that   B,  (t,  ,r,)  -  b^   and   B2(t,  ,r~)  -  b2/ 
because  the  objective  function  which  is  to  be  minimized  is  a 
non-decreasing  function  of   r-^   and   r~   as  are  the  values 
of  the  constraint  expression.   This  indicates  that  an  opti- 
mum solution  can  be  arrived  at  by  varying   r   and   Q   alone 
and  allowing   r,   and   r2   to  assume  the  smallest  integer 
values  such  that  the  constraint  expressions  are  met. 
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It  is  reasonable  to  suppose  that  due  to  the  nature  of 
the  computations  involved  in  the  solution  of  the  program, 
for  example,  summing  various  expressions  which  are  easily 
located  in  tables,  a  high  speed,  digital  computer  would  be 
employed.   While  knowledge  of  convexity  can  be  extremely 
helpful  in  locating  optimal  solutions  and  in  fact  assures 
that  any  local  optimum  is  global  for  the  problem,  the 
limited  nature  of  the  search  region  and  certain  special 
features  of  the  program  give  grounds  for  optimism  that  the 
program  can  be  solved  rapidly. 

It  should  be  noted  that  the  bulk  of  the  computations 
involved  in  a  single  iteration  of  a  computer  search  routine 
involve  determination  of  the   t.  s   as  functions  of   r   and 
Q.   It  is  probably  the  case  that  many  of  the  terms  in  the 
expression  for   tj_   could  be  shown  to  be  so  small  relative 
to  the  others  that  they  could  be  disregarded  without  loss  of 
essential  accuracy.   Central  to  the  computation  of   tj   is 
the  expression 

P(S^s)  =   V    P(S^s  [  H=-h)   P(H=-h). 
h^O 

While  at  first  glance  this  appears  to  involve  formidable 
computations,  the  physical  interpretation  of   H   should  not 
be  overlooked.   It  is  highly  doubtful  that  any  optimum  solu- 
tion to  the  program  would  be  characterized  by  values  of   r 
and   Q  which  would  produce  a  significant  probability  that 
H  be  smaller  than,  for  example,  -3  or  -4,  depending  on  the 
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program  parameters  involved.   Certainly  large  probabilities 
of  small  values  of   H  would  increase   t^   to  a  point  where 
the  values  of  r±      and   r2   required  to  fulfill  the  con- 
straints would  be  prohibitively  large,  especially  when  as  is 
usually  the  case  the  holding  costs  at  the  bases  are  larger 
than  those  at  depot.   In  this  example  it  follows  that  compu- 
tation of   P(S-s)   to  a  degree  of  accuracy  sufficient  for 
the  purposes  of  the  program  would  involve  summing  over  a 
very  few  values  of  the  dummy  variable   h   since  the  remain- 
ing terms  would  be  extremely  small. 

Consider  also  the  regions  of   r   and   Q  which  should  be 
searched.   There  is  reason  to  suppose  that  the  ranges  of 
values  of   r   and   Q  which  produce  minimum  expected  cost 
solutions  to  the  program  are  restricted.   Certainly  for  some 
value  of   r  (even  if  Q=l)  the  contribution  made  by  adding  an 
additional  item  to  the  reorder  point  does  not  significantly 
decrease  the  values  of  the  .  t^'s   and,  therefore,  larger 
values  of   r   need  not  be  considered.   A  similar  argument 
can  be  made  for   Q.   There  is  a  value  of   Q   such  that  the 
benefit  in  decreased  ordering  cost  gained  by  increasing  the 
value  of   Q   does  not  offset  the  penalty  incurred  by 
increased  holding  costs  which  also  result.   These  border 
values  on  the  regions  which  are  searched  for  optimum  solu- 
tions to  the  program  are,  of  course,  extremely  sensitive 
to  the  program  parameter  values. 
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V.   APPLICATIONS  OF  THE  MODEL 

As  noted  when  developing  the  model  it  is  often  the 
objective  of  a  constrained  optimization  problem  to  develop  a 
group  of  alternatives  from  which  a  decision  maker  can  select 
a  "best"  element  which  reflects  his  feelings  about  the  trade 
off  between  comparable  and  non -comparable  costs.   In  this 
thesis  an  analysis  has  produced  a  program  which  provides  a 
method  to  obtain  points  in  a  space  of  costs  and  measures  of 
effectiveness.   For  a  particular  inventory  item  and  situa- 
tion such  a  set  of  points  might  be  as  shown  in  Fig.  2  where 
bj_   represents  the  probability  that   Base  i  is  able  to  fill 
an  order  on  demand,  and  associated  with  each  set  of  values 
of  b-^   and  b2   is  a  cost  rounded  to  the  nearest  thousand 
dollars. 
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Figure  2.   A  Representative  Set  of  Decision  Alternatives 
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While  it  is  not  the  goal  of  this  thesis  to  investigate  the 
various  methods  which  can  be  used  to  develop  a  cost- 
effectiveness  surface,  the  following  points  are  worthy  of 
note: 

The  logic  of  the  physical  realities  of  the  model 
dictates  that  costs  will  be  non-decreasing  in  bj_   and  b2- 
If  this  were  not  the  case  it  would  be  possible  to  "buy"  a 
greater  amount  of  effectiveness  for  less  than  is  needed  to 
pay  for  some  lesser  degree  of  protection  which  is  contradic- 
tory to  the  minimum  cost  result  of  the  efficient  solution. 

In  many  situations  it  will  only  be  necessary  to  develop 
minimum  cost  solutions  for  small  ranges  of  possible  values 
of  b-|_   and   b2«   The  reason  for  this  is  that  in  many 
instances  the  decision  maker  can  specify  in  advance,  without 
immediate  knowledge  of  the  consequences  in  terms  of  effec- 
tiveness, the  maximum  expected  cost  which  can  be  incurred 
and/or  the  general  range  of  effectiveness  he  wishes  to 
achieve.   It  may  be  noted  that  certain  ranges  of  effective- 
ness might  be  of  particular  interest  and  that  for  those 
areas  the  effect  of  a  smaller  change  in  the  values  of  b^ 
and  lo>2      could  be  sought;  therefore,  it  would  be  necessary 
to  search  along  a  finer  grid  of  values  for  b^   and  b2   at 
least  for  some  particular  segment  of  their  ranges.   It  can 
be  said,  therefore,  that  the  ranges  of  b-.   and  b2   to  be 
searched  can  be  narrowed  considerably  with  knowledge  of  the 
decision  maker's  general  policy  and  the  physical  realities 
of  the  system  to  which  the  model  is  applied. 
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It  might  have  been  noted  that  terms  in  the  objective 
equation  involving  the  repair  facilities'  holding  costs 
were  independent  of  the  program  variables  and,  therefore, 
could  have  been  deleted.   These  terms  were  included,  how- 
ever, in  order  to  show  the  ease  with  which  this  model  and 
program  can  be  employed  to  discover  not  only  efficient 
inventory  ordering  policies  but  also  efficient  inventory 
ordering  and  repair  facility  management  policies.   Suppose, 
for  example,  it  is  possible,  by  incurring  an  extra  cost,  to 
decrease  the  repair  time  at  the  depot  by  some  fraction  of 
the  original  repair  time,  that  is  there  exists  a  function 
C(rt)  defined  for   rt   in  some  interval  j a,bj  ,  which  is 
decreasing  in   rt.   It  is  then  possible  to  specify  a  new 
program: 

Minimize        K(rx.  r2,  r,  Q,  rt)  +  c  (rt) 

subject  to      Bi(ti,  r±)  ^  h± 

rt   >   a 
-rt   £-    -b 
i   =•  1,  2 

where   r^,  r2,  r   and   Q   are  non- negative 
integers  and    rt   is  some  real  number. 

Similar  cost  terms  and  constraints  can  be  added  for  base 
repairs  or  any  combination  of  repairs  where  the  repair  times 
can  be  considered  variables.   In  the  case  of  the  bases, 
however,  the  additional  costs  need  only  be  considered  to 
vary  with  the  expected  base  repair  times.   It  is  then  pos- 
sible as  above  to  present:  the  decision  maker  with  a  revised 
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set  of  costs  corresponding  to  the  various  degrees  of  effec- 
tiveness achieved. 

It  may  also  be  noted  that  while  the  model  configuration 
developed  has  two  bases,  there  is  no  reason  why  systems  with 
any  number  of  lower  echelon  units  cannot  be  modeled.   The 
only  changes  required  in  the  program  formulation  would  be 
to  increase  the  number  of  cost  terms  in  the  objective 
function  and  the  number  of  constraints  by  a  number  of  terms 
equal  to  the  additional  bases  modeled.    it  may  be  observed 
that  the  model  specified  can  be  used  to  describe  a  multi- 
echelon  system  where  item  repair  is  not  considered.   This  is 
possible  by  stipulating  the  probabilities  of  repair  at  the 
various  facilities  to  be  zero. 
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